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We used Geometric Algebra to compute the paths of skew rays in a cylindrical, step-index mul¬ 
timode optical fiber. To do this, we used the vector addition form for the law of propagation, the 
exponential of an imaginary vector form for the law of refraction, and the juxtaposed vector product 
form for the law of reflection. In particular, the exponential forms of the vector rotations enables us 
to take advantage of the addition or subtraction of exponential arguments of two rotated vectors in 
the derivation of the ray tracing invariants in cylindrical and spherical coordinates. We showed that 
the light rays inside the optical fiber trace a polygonal helical path characterized by three invariants 
that relate successive reflections inside the fiber: the ray path distance, the difference in axial dis¬ 
tances, and the difference in the azimuthal angles. We also rederived the known generalized formula 
for the numerical aperture for skew rays, which simplifies to the standard form for meridional rays. 


I. INTRODUCTION 

Optical fibers are waveguides used for optical commu¬ 
nication as first illustrated by Tyndall in 1870[T]. Single¬ 
mode fibers have radii of about 8 microns, while multi- 
mode fibers have about 60 microns [5]. The description 
of light propagation in single-mode fibers require wave 
optics, while that in multimode fibers, the ray approxi¬ 
mation is sufficient. [3] In this paper, we shall only talk 
of multimode fibers using ray or geometric optics, i.e, no 
diffraction, interference, or polarization. 

Optical fibers can be classified depending on the func¬ 
tional dependence of the refractive index on the fiber 
radius: step-index or graded-index (GRIN). Step-index 
can be multi-step [4] while, graded-index fibers can have 
refractive index functions that are parabolic. Here, we 
shall focus only on step-index fibers, where the refractive 
index of the core is constant and light rays are guided 
through total internal reflection. 

If the rays lie only on a single plane, the rays are said to 
be meridional; if the rays trace a discrete helix, the rays 
are skew. Most textbooks describe only the meridional 
rays and mention skew rays only in passing. The reason 
for this is that the mathematics used to describe skew 
rays is difficult. m 

The study of geometric optics ray tracing for skew rays 
is still an important problem for two reasons. The first 
is pedagogical: skew ray tracing for laser light in large 
fibers, e.g., HeNe laser in 1 cm diameter glass fiber, can 
be used to verify the helical properties of skew rays in a 
student laboratory. The second is practical: skew rays 
tracing equations can serve as the starting point for the 
analysis of electric fields at the fiber walls [5] and can 
provide the basis of ray tracing algorithms for optical 
fibers |9], which may even include the Goos-Hanchen ef¬ 
fect for reflection [TU]. 
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To construct a ray tracing algorithm, we need to de¬ 
scribe rays mathematically. We can do this in three ways. 

The first way is by labeling points, lines, and angles, 
then use geometry and trigonometry to determine the 
scalar invariants, e.g., the axial path length between re¬ 
flections, axial angle of incident and reflected rays, and 
the numerical aperture. But the system for naming these 
geometrical quantities vary from author to author, mak¬ 
ing it difficult to compare results. EnnHH] To solve the 
problem of the angle naming conventions, we define the 
directions of the vectors in terms of the polar angles 9 
and the azimuthal angle (j) in spherical coordinates, dis¬ 
tinguished only by their subscripts, e.g. 9cr and for 
the incident ray cr and for the normal vector rj. 

The second way is to use vectors for ray tracing. In 
general, we need five vectors: the initial position of the 
ray r, the final position of the ray r', the propagation 
vector cr, the reflected vector cr', and the vector normal 
to the surface rj'. The task of ray tracing then is to 
relate these variables using a set of equations, such as 
those given in Klein and Furtak[T5]. though in a slightly 
different form involving the concavity function Car/k = 
± 1 : 

Tfe-i-i = ffc + Sfe-i-i (la) 

cr^ — Cr/2 ‘^t^arjk'Hk cos/?/j;, (^t)) 

tlk+l^k+1 — tlk^k 

“f Ccrr/A: COSCOS(Ic) 

which corresponds to propagation, reflection, and refrac¬ 
tion. But these ray tracing equations are difficult to ap¬ 
ply if we wish to compute for analytic solutions, such as 
reflected ray and final position after a certain number of 
reflections within the fiber. That is why ray tracing algo¬ 
rithms do not use analytical solutions but rather iterative 
computations using the equations described |!I1 IT(T|. 

The third way is to use vectors in ray tracing but 
within the framework of the Eikonal method. [ig Hz] 
Here, the position R of the ray satisfies a vector differ- 
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ential equation: 


II. GEOMETRIC ALGEBRA FOR GEOMETRIC 
OPTICS 


d 

ds 



= Vn, 


( 2 ) 


A. Geometric Algebra 


where s is the propagation distance and n is the refractive 
index. This method lends itself readily to the analysis of 
invariants in step-index fibers as shown by Zubia et al. 
[3]. This method may also be applied to graded-index 
fibers using computational methods O [HI [19]. This dif¬ 
ferential equation approach may be powerful and lends 
itself well to computational methods, but it does not ex¬ 
ploit the simple vector geometry of rays for step-index 
fibers. 


In this paper, we shall not use the Eikonal method 
nor the reflection and refraction equations in Eqs. (lb) 
and (Ic). Instead, we shall construct the ray tracing 
equations using the exponential rotational operators and 
direct vector products in Geometric Algebra: 


ab = a b-|-iax b, (3a) 

g*ne _ 0 + in sin 9, (3b) 


which are the Pauli identity[5ni [2T] and Euler’s Theo¬ 
rem, respectively. For example, the law of reflection and 
refraction can be expressed as 


^fc+1 — 'nk^k'Hk-} 

(4a) 

(Tk+l — + 

(4b) 


The direct vector product form of the law of 
reflection in Eq- (4a) and the exponential rota¬ 

tion form of the law of refraction [Mj in Eq. (4b I were al¬ 
ready known before. These laws and their variants were 
used in geometric optics for the derivation of ray tracing 
equations for spherical lenses and mirrors for different 
ray cases: finite skew, paraxial skew, finite meridional, 
and paraxial meridional [251 [26] . Equations (4a) and (4b) 


were also used before in the derivation of skew ray tracing 
equations in optical fibers [T^: however, the treatment 
in this three-page extended conference abstract is lim¬ 
ited and does not include the discussion of the numerical 
aperture. We shall correct these in this paper. 


We shall divide the paper into six sections. Section|T|is 
Introduction. In Section jll] we shall describe the basics 
of Geometric Algebra and how it can be used for vector 
rotations in different coordinate systems. We shall use 
this algebra to describe the Laws of Propagation, Reflec¬ 
tion, and Refraction in Geometric Optics. In Section jllTl 
we shall derive the distances and angles of propagation 
of the light rays inside the fiber and compute the fiber’s 
numerical aperture. In Section |IV[ w e shall summarize 
the ray tracing invariants. Section is Conclusions. 


1. Vectors and Imaginary Numbers 

Geometric Algebra known as Pauli Algebra. [20] 
is generated by three spatial unit vectors ei, 62 , and 63 
which satisfy the orthonormality axioms m- 

= 1 , (5a) 

Bjek = -ekBj. (5b) 

That is, the square of a vector simply is unity and the 
juxtaposition multiplication product of two distinct unit 
vectors anti-commute with each other. 

We can generalize these multiplication rules for two ar¬ 
bitrary vectors a and b in 3D using the Pauli Identity [20l 

mi 


ab = a • b -I- i{a x b), (6) 

where z = 616263 . It can be shown that i is an imaginary 
number that commutes with both scalars and vectors. 

If we express a as the sum of its perpendicular and 
parallel vector components with respect to b, then 

a = aj_-kaj|. (7) 

Using the Pauli Identity in Eq. (|^, we can show that 

ayb = ay • b = bay, ( 8 a) 

axb = ax X b = —bax. ( 8 b) 

Equation ([^ is the 3D generalization of the orthonormal¬ 
ity axiom in Eq. ([^. 


2. Rotations 


Let 0 be a scalar and n be a unit vector. Since (zn)^ = 
— 1, then we may use Euler’s Theorem: 

g*ns _ COS0 -I- znsin0. (9) 


Now, suppose ax and ay are the components of a per¬ 
pendicular and parallel to n. Multiplying these from 
the left of Eq. ([^, and using the commutation and anti¬ 
commutation rules in Eq. ([^, we obtain 


anU"® = U"«an 


axe™® =e-™®ax. 


(10a) 

(10b) 


Notice that a change in sign of the exponential’s argu¬ 
ment occurs only when the exponential is multiplied to 


ax. 


If we expand the left side of Eq. (10b), we get 


ax e™® = ax (cos d -I- zn sin d). 


( 11 ) 
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FIG. 1: Relations of the Cylindrical Coordinate System 
as compared with Cartesian and Spherical Polar. 


FIG. 2: A vector Sk+i^k+i is defined using two vectors 
with the same origin point, and r^+i. 


Distributing the terms, and using the Pauli identity in 
Eq. (|^, we obtain 


a_Le*"^ = ax cos 9 — (ax x n) sinf 


( 12 ) 


since ax • n = 0. Geometrically, axe*"® is the vector ax 
rotated counterclockwise about a unit vector n. 


3. Coordinate Systems 

We claim that the unit radial vector e,. in spherical 
coordinates is given by [53] 


63 e“2egie3<)!./2^ 


(13) 


To show this, we note that from Eq. (12), 


636 *®^® _ ggg 6 * + 61 sin 9, 


so that Eq. (13) may be written as 

(63 cos 9 + 61 sin0) 




(15) 


Distributing the terms and using the relations in Eq. (10 1 , 
we get 


= 63 cos 9 + Oi sin 9 


(16) 


Hence, 


= 61 sin 9 cos ^ + 62 sin ^ sin 0 + 63 cos 9. (17) 


Equation ( |17| ) is the representation of the unit spherical 
radial vector in rectangular coordinates, as shown in 
Figure 12 


B. Geometric Optics 

1. Law of Propagation 

If a light ray moves from its initial position at the 
interface by a distance Sfc+i in the direction of the 
unit vector <Tk+i, then the ray’s final position r^+i at the 
{k + 1)*^ interface is [24] 


rfc+l — Tfc + Sk+l CTk+l- 


(18) 


(see Fig. [^. Note that we can express r^xi, as 
the sum of 63 and eie*®^*^ as given in Eq. (16) : 


Tfe = Xfces + r^;Sin 6 <feele*®’*‘^^ (19a) 

rfc+i = Zfc+iea + r^+i sin6<fe+ieie*®’5'^''+b (19b) 

(Tk+i = cos 6 »<^(fc+i)e 3 +sin 6 ',^(fe+i)eie“’^'^''''+b (19c) 


(14) then Eq. (18) can be separated into simultaneous equa¬ 


tions for each component. 

Cfe+i sin 9k+i cos (j)k+i = sin 9^ cos cj)k+ 

afc+i sin ^cr(fc-t-i) cos (/>cr(fc+i) 7 

( 20 a) 

Tk+i sin 9k+i sin (/)/c+i = sin 9k sin (j)k+ 

'Sfc-t-l Slci ^cr(fc-t-l) 4^(T{k+l) 7 

( 20 b) 

cfe+i cos 9k+i = rk cos 9k + Sk+i cos 9cr(k+i) ■ 

( 20 c) 

2. Law of Reflection 

Let (Tk be the direction of the ray as it hits the 
interface with a unit normal vector rjk- The new direction 
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(Tk+i after it hits the interface is 


(Tfe+i = 


( 21 ) 


That is, the reflected ray is the negative of the 

incident ray cr^ rotated counterclockwise about rik by 
an angle tt. Expanding the exponentials using the Euler 
identity Eq. (§, then Eq. ( [^ reduces to 




( 22 ) 


which is a known form for the reflection law in Geometric 
Algebra. [3S] 


If we rewrite the product rj^cr^ in Eq. (22 1 using the 
Pauli Identity in Eq. ^ and distribute the rightmost r]k, 
we get 


cTk+i = -{Vk ■ crk)rik + iiVk x (Tk)rik. 


(23) 


where cr^ is the incident vector, (Tk+i is the refracted 
vector, Ca-rik is the concavity function in Eq. (27), and 
e„r]k is the axis of rotation 


-'(7r)k 


cTfc X r}k 

Wk X rjfcl 


(31) 


with rjk is the unit vector normal to the surface. Geomet¬ 
rically, Eq. (30) says that the refracted ray is the incident 


ray rotated either clockwise (if the concavity function is 
negative/interface is convex) or counterclockwise (if the 
concavity function is positive/interface is concave) about 
^arjk by the difference of the incident angle and the re¬ 
fracted angle, /3 —/?'. We will use primed variables for the 
refracted values. Also, we shall mostly deal with concave 
surfaces inside the fiber, so that we can immediately set 

^CTTjk = “t“l. 


Applying the Pauli identity again for the second term 
results to 


cTfe+i = -r]k{r]k ■ (Tk) + [r]k x crk) x qk, (24) 


since i{qk x crk) ■ Vk = 0 . Einally, expanding the triple 
cross product on right hand side, Eq. (24) becomes 


(Tk+i =-^VkiVk ■ (^k) + f^k- (25) 


Let us define Pk as the angle of incidence, so that 


Vk ' — ^ar]k COS /I/j;, (26) 


where Ca-rik is the concavity function. 


Ccrr/fc 


^k ' 'Hk 
k/c • r7fc| ’ 


whose values are ±1. Thus, Eq. (25) reduces to 


crk-\-l — ‘^Vk cos I3k “t" CTk , 


(27) 


(28) 


which is the same expression for the law of reflection in 
Klein and Furtak. [151 IMl 


3. Law of Refraction 
The Law of Refraction is given by 


III. SKEW RAYS IN OPTICAL FIBERS 

A. Refraction at Point Pi 

Let (Ti be the direction of the initial ray as it strikes 
the flat end of the fiber: 

(Ti = sin^CTi -|- cos9^. 63 , (32) 

which is similar in form to Eq. ( |19c[ ) . If the ray (Ti strikes 
the fiber at the position r^, then the outward normal 
vector at the surface is 


I?* = -ea- 


(33) 


Multiplying Eq. (32) and Eq. (33), and using the ex¬ 


ponential identities in Eq. (1^, we get 


CiVi = sin 9ai *626 


?,©3 (j^CT 


— cos . 


(34) 


Expanding the left-hand side of Eq. (34) using the Pauli 
Identity in Eq. (§, and separating the scalar and vector 
parts, we obtain 


Ci ■ Vi = ~ COS 9ai, (35a) 

(Ti X qi = sin0(T, 62 (35b) 


Substituting these into the definition of the the concavity 
function and the rotation axis in Eqs. (27) and (31), we 
arrive at 


n sin ff = n' sin j3 ', 


(29) 


where n and n' are the indices of refraction at both sides 
of the interface, while /3 and P' are the angles of incidence 
and refraction. 

The refraction law in Eq. (29) may also be reformulated 
in geometric algebra as follows 


cos 9^. 

= "|COS 0 <,. I = 
ear,i = 626 *®"'^-'-, 


(36a) 

(36b) 


since 9^. > 0 . 

We can find the angle Pi between the unit vectors (Ti 
and qi using their dot product: 


^fc +1 — ^k 


^iCa-qGa-q{l3 j3 ) 


(30) 


CTi ■ Vi = COS Pi = cos 9ai , (37) 
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(a) The incident ray ai strikes the point Pi entrance end of the 
fiber where the normal vector is rji. The angle of incidence is 
Pi and the angle of refraction is PP From point Pi, the light 
ray enters the fiber and strikes the cylindrical walls at points 
Po, Pi, and P 2 . 


Pi 



(b) Top view of Fig. (3a|, with the points Pg, Pi, and Pi as the 
projections of the points Pg, Pi, and P 2 on the fiber’s entrance. 


FIG. 3: The propagation of the light ray inside the 
optical fiber in perspective view and top view 




FIG. 4: Vector Tq can be expressed as the sum of a 
rotating vector on the xy plane, Rr]o, and a vertical 
component, ^ 063 . 


so that Eq. (40) reduces to 


Co = 636 


Expanding the exponential in Eq. (421, we get 


(To = sin/3' + cos/3' 63 . 


(42) 

(43) 


On the other hand, using the expansion of (Tq in 
Eq. (19c), we have 


(To = cosS^eo + sin6i_eie“’^‘^‘'o. 


(44) 


Thus, equating Eqs. (43) and (44), we get 


= /3l, 


(45a) 

(45b) 


which is the same relations in once we have changed 
the subscripts to match their notation. 


B. Propagation from Pi to Pg 


so that 

= (38) 

Using Snell’s Law, the refracted angle /3' is given by 

^'= sin"i (^^sin/3i) , (39) 

where we set n = 1 for air (see Fig. |^. 

The law of refraction as given by Eq. ([M| for fc = 0 is 


0-0 = 


(40) 


where Co-,,! = 1. We can express cr^ as a rotation of 63 
about the axis Sa-rii by a clockwise angle Pi, 


We define the entry point of the light ray when it in¬ 
tersects the base of the cylinder as r^: 

Xi = Pj = Pi cos piei + sin 0 ^ 62 . (46) 

From position r^, the ray moves in the direction of ctq de¬ 
fined in Eq. (44), until it strikes the walls of the cylinder 


at position rg, so that 

rg = Sgcrg -|- r^. 


(47) 


where sg is the distance of propagation. If we define the 
radius of the cylinder to be R, then the position rg can 
be written as 


(Ti — 636 


(41) 


rg = -b zgeg. 


(48) 
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where is the azimuthal angle and zq is the distance 
along the axis, as illustrated in Figure Substituting 
Eq. (47|) into Eq. Esl), we get 


soCTo + Vi = + zges. 


(49) 


Notice that there are three unknowns in this equation: 
So, and zq- 

Let us separate the axial and radial components of 
Eq. (49) to obtain 


So sm^o-geie 


socos^o-oeg = zges, 
^0 + Pi eie" 


_1_ n 


(50a) 

(50b) 


Equation (50a I gives us our first equation relating sq and 
Zq: 



FIG. 5: If we extend so+crg backwards, we obtain a 
vector of length so+, which will terminate at the fiber’s 
wall, extended in the —eo direction. 


zo = So cos 00 


(51) 


Now, squaring both sides of Eq. (50b) and using the 
exponential identities in Eq. (10), we get 

Sg sin^ 0^0 +Pi +2so/9i sin0„Q cos((f>i - . (52) 

Solving for sq sin O^jq using the quadratic formula results 
to 


So sin0^g = ±^B? - pI - ^^o) 

- Pi cos{(j)i -(l>ao), 


(53) 


after simplifying the terms. 

Notice that there are two possible solutions: 


1 


so± = 


sm( 


0-0 


± 


\Jr‘^-p1 sin^(((ii - 0 oo) 


- Pi COs{(l)i - (jica) 


(54a) 


The distance so+ is the distance traveled by the light ray 
inside the fiber from the fiber’s entrance to the cylindrical 
interface, while the distance So_ is the distance traveled 
by the light ray if it propagates opposite to the direction 
(To. The difference between these two distances is 


that for our case, so- is a virtual backward propagation 
from the fiber entrance. 


Lastly, we solve for (j)rjg. Separating Eq. (50b) into the 
components of ei and e 2 , we get 


So sin 0cr„ cos (pao + Pi cos 4>i = R cos , 
So sin 00-0 sin + p^ cos pi = R sin p^,^. 


Dividing Eq. (58a) by Eq. (58b), we obtain 

Uvo 


tan pi^Q — 


(58a) 

(58b) 


(59) 


^Vo 


where 


Xr)o — 


So sin cos Pag + Pi COS pi, (60a) 
y-qo = So Sin0<^^ s'mpao + Pi sin Pi- 


(60b) 


We note that So sin 0^0 is given in Eq. (53), so that the 


rectangular coordinates ana 
rjo are already defined in terms of the initial parameters. 
Thus, the cylindrical coordinates of the normal vector rjQ 
can now be determined, since the radius R of the fiber is 


constant and the azimuthal angle prjg is given by Eq. (591. 


C. Refraction at Pqi Numerical Aperture 


so+ - So- = . ^ J- Pi - <(>cTo), (55) 

Sin c7^q 

which is a new result (see Fig. [^. If pi = 0, then we get 
the meridional case: 


so+ — Sq- 


2R 

sin 9cro 


(56) 


At point Pq, on the fiber’s cylindrical wall, the incident 
ray is (Tq defined in Eq. (44), and the normal vector is 


r]o = 


(61) 


Their product is 

cToVo = sin 0 croe*®^^‘^''o“'^''o) — i cos 0^0 026 *®’^, (62) 


On the other hand, if pi = R, 


So+ ~ So- 


2i?| cos{pi — 
sin 0crg 



(57) 


which is similar to the form of the wall-to-wall propaga¬ 
tion distance given by Cozannet and Treheux|28j. except 


where we used the exponential indentities in Eq. ( |10[ ). 
Separating the scalar and imaginary vector parts, we get 

(To • J7o = sin 0,jg cos{pag - Pqo)i (63) 

o-Q X m = ei cos 0cro sin pn^ - 02 cos 0o-o cos pr,o 

+ 03 sin 0^,0 sin((;icro - (Pvo)^ (64) 
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after removing the imaginary number i in the second 
equation. If we define ip' as the angle between erg and 
?7o, then Eq. (IMl) yields 


cosip' = sin6>„^ cos{(pao - (pno)^ 


(65) 


which is essentially the same as that of Potter [TT] and 
Senior 

Since ip' is also the angle of incidence at point Pg, then 
by Snell’s Law we have 


n' sinip' = nsin'0, 


( 66 ) 


where ip is the angle of refraction outside at point Pg. 
For the ray to be trapped, we set i/i = 7r/2, so that 


sm'Ip = 


Combining Eq. (65) and Eq. (67), 


(67) 


1 = ^ + sin^ cos^{cp„^ -(pr,o), 


( 68 ) 


then solving for sin^o-o, we get 


sin 0 ^„ = 


\/ n '2 — 


z'|cos(()),^„ -(pno)\' 


(69) 


Now, the numerical aperture NA of an optical fiber is 
defined as 


NA = n' sin 0a 


We combine Eqs. (69) and (70) to get 


NA = 


\/ n '2 — rfi 
COs((^ao 


(70) 


(71) 


Equation (71) is the expression for the numerical aper¬ 
ture for skew rays in optical fibers as given by Potter pdl 
[22] and Senior [T^. 

For the meridional approximation, we set (prjo = (pao, 
so that Eq. 0 reduces to the known standard form for 
the numerical aperture as given by Klein and Furtak[T5] 
and PotterflT]: 


NA = \/n'2 -nA (72) 

If we expand the denominator in Eq. we get 
cos{(p„o - (pT,o) = cos pao COS pr/o + sinpag sinp^jg. (73) 
Since 


cosprjg = Xr,o/R, (74a) 

sinpno =yr,o/R, (74b) 


then using Eq. (60) together with trigonometric indenti- 
ties, we obtain 


cos{pao 


<^^, 0 ) = ^ sin Oao + ^ cos{pi -pao)- (75) 



FIG. 6: The skew factor x as a function of the ratios of 
the relative radii p/R and the azimuthal angular 
difference pi — pa^ 


Substituting Eq. (75) back into Eq. 0 results to 

y/n '2 — 71,2 


NA = 


^ sin 9ao + ^ cos{pi - pao) 


(76) 


Combining this with the expression for sg in Eq. (54a) 
and using the identity pao = Pa-i from Eq. (45b), we 
arrive at 


NA = xV'n-'^ ~ 'cR': 


(77) 


where 


X = 


Pi 


(78) 


1 - 


(Pai) 


Equation (77) forms the expression for the numerical 
aperture NW^ for the skew rays in optical fibers. Here 
we introduce the skew factor y which is defined in terms 
of the initial ray parameters at the fiber entrance: the po¬ 
sition Ti = [pi,pi,Q) in cylindrical coordinates where the 
incident ray er^ = (l,pa^,0a^) in spherical coordinates. 
Notice that the skew factor depends on the azimuthal 
angle pa- of the incident ray, but not on the polar angle 

The expression for the numerical aperture NA in terms 
of the skew factor % is a new result. To arrive at the 
meridional case in Eq. (72), we set the radial distance 


Pi = 0, to get the skew factor x = 1- On the other hand, 
if we set the initial position at the edge of the fiber, i.e. 
Pi = R, we get 


1 


X = 


y/cOS^{pao - Pvo)' 


(79) 


which leads to Eq. 0, the numerical aperture for skew 
rays given in the literature. Fig. shows the plot of 
the skew factor x as a function of the normalized radial 
distance Pi/R and the difference of the azimuthal angles 





































D. Reflection at Pr 


0 


At position Tq on the interface, the light ray’s propa¬ 
gation vector is ctq and the normal vector to the interface 
is r/o: 


cTo = cosOa-oes -I-sinSo-oeie^’^'^'^o, 


rjo = 


(80a) 

(80b) 


To obtain the new propagation vector cri after reflection, 
we use the law of reflection in Eq. (22): 


fi = —'noO’oVo- 


Substit uting the expressions for t/q in Eq. (80b), and ctq 
in Eq. (80a I into Eq. (81), yields 


(Ti = — Sin 


+cos6»„ 


^3 5 


(82) 


after distributing the terms and simplifying the result. 
Eq. (82) is the desired equation for the new propagation 


vector (Ti. 


Using the definition of cri as written in Eq. (19c) 
(Ti = sin 0(71 -|- cos den 63 , 


(83) 


and comparing this with the expression for cri in Eq. (83), 
we obtain El 


^CTo ; 


— 4 ^cr 0 - 


(84a) 

(84b) 


Notice that Eq. (84a) is not proper as we previously de¬ 


fined 9 to have values between 0 to tt. 

In order to avoid having to write the negative sign 
before the 9, we can use the following relation derived 
from the Euler identity: 


g2e3(±7i-) _ Qos(± 7 r) -I- ie 3 sin(± 7 r) = —1, 


so that Eq. (82) becomes 

cTi = sin0,7g 

Hence, 




0±7r) 


-b COS 0(70 eg. 


(85) 

( 86 ) 


(87a) 

(87b) 


4^(70 i 

which is similar to what we have derived in a previous 


paper[T3]. Notice that Eq. (87a) is now in the proper 
form of 0, except that the sign of tt in Eq. (87b) must be 
determined for individual cases. 


E. Propagation from Pq to Pi 


From the propagation law in Eq. (18), we know that 


the final position of the ray depends on its initial position 

ro: 

( 88 ) 


where si and cri is the ray’s propagating distance and 
direction, respectively. This is shown in Figure If we 
use the expression of Tq in Eq. (|48|) and express rg in a 


similar form, then Eq. ( 88 ) becomes 

siCTi -b -b ^063 = -b 2 : 163 . (89) 


Separating Eq. (89) into its radial and axial components, 
we get 

Si sin 0 ( 7 ^ -b i?eie*®^'^’>o = i?eie*®^'^’>i, (90a) 

Si cos 0(71 63 -b zoes = zie 3 . (90b) 


(81) Notice that there are three unknowns in Eq. (90): Si, 


0r)i J z:i. 

To solve for the propagation distance si, we first square 
both sides of Eq. (|90a|): 


Si sin^ 0(71 -bi?^ -b 2 sii?sin 0 ( 7 j cos{(j)no — 4>ai) = ■ (91) 

Hence, 

2Rcos{(j)j^g - (()(7i) 


Si = — 


sin 0(7 


(92) 


To solve for (f>rin we expand Eq. (90a) into its ei and 
62 components: 

R cos ((> 7,1 = Si sin 9^^ cos 4>ai + R cos (pr^o, (93a) 

RsiTKprfi = Si sin0(71 sin(()(7i -b Rshre^r^g. (93b) 

Dividing the two equations and using the expression for 
Si in Eq. (92), we get 


_ sin cos 2 (/)( 7 i - cos ((iim sin 

hdblKpjj, — - - -——-;- - -; 

— COS COS 2(j)(j^ — sin sin 20^^ 
after using the double angle identities. Hence, 

-sin(2((i,7i - (j>r,o) 


tanc^ij^i — 

- cos(2 (()(7 j - . 


'Vo 


y 


(94) 


(95) 


Note that we retained the negative signs in both the nu¬ 
merator and denominator so that we can properly deter¬ 
mine the angle 0 ,;i which is from 0 to 27r. 


IV. SKEW RAY INVARIANTS 

A. Propagation Distance Between Reflections 

The procedure for computing the propagation of the 
light ray from Pi to P 2 is the same as the procedure 
from Pq to Pi. Therefore, we can simply change the 
subscripts to obtain the following relations for the prop¬ 
agation distance S 2 , the propagation polar angle 9 „^, and 
the propagation azimuthal angle (j)a 2 '- 


S2 = 


9(71 — 


2Pcos((^,)i - ^(73) 

sin 9rro 


<71 , 

Ul 


0cr2 — ^^771 4^(71 2b TT. 


(96a) 

(96b) 

(96c) 


n = siCTi -b ro 
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Taking the difference of Eq. (102 1 and Eq. (103), and 
imposing S 2 = si in Eq. (100) and 9^2 = 9a^ in Eq. (96b), 
then 


l\z = Z 2 - Zi = Zy- Zo¬ 


rn) 


Thus, the axial distance Az between wall-to-wall propa¬ 
gations, defined as the pitch, is constant. 

If we substitute the expression for 21 in Eq. (102) and 
Si in Eq. (92), we get 


Az = -2R cos((/>no - ) cot , 


(105) 


which is similar in form to that given in Eq. (22) of 
Cozannet and Treheux|28j. Alternatively, we can express 
Az in terms of si in Eq. (92) to get 


Az = Si cos 9„ 


(106) 


FIG. 7: The light ray traveled from point rg to point ri 
on the fiber’s walls, at a distance si, along the direction 
(Ti. The normal vectors at the two points are rio and rji. 


which is an invariant given in Love and Snyder[3]. Note 
that this invariant reduces only to the ray half-period if 
the rays are meridional. 


Notice that the angle 6^ remains constant through reflec¬ 
tions. Our aim is to show that the propagation distance 
s is also constant between reflections, i.e., si = S 2 . 


If we substitute Eqs. (96b) and (96c) into Eq. (96a), 
we get 


S2 = 


2R cos{-(j)j^^ -k (t)aj 
sin 9 c,, ’ 


(97) 


since cos((/) — tt) = —cos(p. Expanding Eq. (97) using a 
trigonometric identity, and using the expression for the 
sines and cosines of in Eq. (95), 


sm(j)r,2 = -sin(2(;icri - (l)r,o), 
cos — COs(2^(j^ ), 


we obtain 


S 2 = 


2 flcos((()^P - 
sin( 6 *^J 


(98) 

(99) 


( 100 ) 


Comparing this with the expression for si in Eq. (92), 
we arrive at 


Si = S 2 . 


( 101 ) 


Thus, the propagation distance between two reflections 
is constant. 


B. Axial Distance Between Reflections 


The axial position zi of point Pi is given in Eq. (90b): 

zi = Si cos^cri + ^0- (102) 


C. Change in Azimuthal Angle Between 
Reflections 


We know from Eq. (95) that the azimuthal angle (pri^ 
at point Pi is given by 


so that 


V-cos( 2 (/)<,^ - 


— 2 ()>(j^ (j^rjQ T kir, 


(107) 


(108) 


where k is - 1 - 1 , — 1 , or 0 . 

To find fc, we substitute the expression for _ in 

Eq. (108), for in Eq. (96c), and sin(0crj in Eq. (96b) 
to Eq. (96a) to obtain 


S2 = 


2i7 cos(2(^fji 0770 207^^ 0(71 T (^ T 1)^) 

sin 9 c,, 


(109) 

Substituting the expression for 07^1 in Eq. (108) again 
and rearranging the terms, we get 

2i?cos(-0(7i + 0770 + (-fcT l)7r) 

S2 = —-----—■ (.IlUj 

sin drr-, 


Comparing this the expression for S 2 in Eq. (100), we 
arrive at = =f 1 , so that Eq. (108) reduces to 




( 111 ) 


Using a similar process, we can show that the form of 
07,2 and 07,3 are similar to 07 ,^, 


Similarly, we can show that for point P 2 , 

Z2 = S2 cos 9cr2 + Zi. (103) 


0772 — 200-2 0 »?i T ( 112 a) 

0^73 “ 200-3 0772 T '^1 ( 112 b) 
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so that 


4‘ri2 4’rii — 2(00-2 (113a) 

0r;3 'Pr]2 2(0o3 0J?2) T (113b) 


Replacing the subscripts of Eq. (96c) to obtain an ex¬ 
pression for 003 ! 


— 20„3 


± tt, 


(114) 


and using the th e expr ession s for 0o3 in Eq. (114) and 
0^2 in Eq- (112a), Eq. (113b) becomes 


— 2(002 'Pvi) T 


(115) 


Comparing Eq. (113aI with Eq. (115), we arrive at 


A0 0^3 0„2 


"V2 


"Vl ■ 


(116) 


Thus, the change in the azimuthal angle A0 between 
consecutive reflections is a constant. 

So now we have two invariants between consecutive 
reflections: the change in azimuthal angle A0 and the 
change in the axial propagation distance Az. From this, 
we conclude that the rays inside a fiber trace a polygonal 
helix. 


V. CONCLUSIONS 

In this paper, we used Geometric Algebra to compute 
the paths of skew rays in a cylindrical, step-index optical 
fiber. To do this, we used the vector addition form for the 
law of propagation, the exponential of an imaginary vec¬ 
tor form for the law of refraction, and the vector product 
form for the law of reflection. In addition, we used the 
spherical angles 9 and 0 to describe the directions of rays 
in space, but expressed in cylindrical coordinates in ex¬ 
ponential form. We showed that the light rays inside the 
optical fiber trace a polygonal helical path characterized 
by three invariants between sucessive reflections: (I) the 


ray path distance, (2) the change in axial distances, and 
(3) the change in the azimuthal angles. We also showed 
that the numerical aperture for skew rays we obtained is 
the same as that of the literature. 

To derive the ray tracing invariants, we did not use 
the reflection and refraction laws in Klein and FurtakflK]. 
Rather, we used two alternative techniques: (a) the Pauli 
Identity which expresses the geometric product of two 
vectors as a sum of their dot and imaginary cross prod¬ 
ucts and (b) the Euler’s Theorem which expresses the 
exponential of an imaginary vector as a sum of the co¬ 
sine of the magnitude of the vector and the product of 
the normalized imaginary vector with the sine of the vec¬ 
tor’s magnitude. These two theorems allow us to express 
vector rotations in exponential form, which enables us to 
take advantage of the addition or subtraction of exponen¬ 
tial arguments of two rotated vectors in the derivation of 
the ray tracing invariants. 

Many of the equations for the invariants were already 
known before, except maybe for the change in the az¬ 
imuthal angle between reflections. Also, we obtained a 
new expression for the numerical aperture, which allows 
the point of entry of light to be an arbitrary point in 
the fiber’s entrance, and not limited to the center or the 
edge as given by previous authors. Finally, our use of 
standard notations for the angles 9 and 0 for cylindri¬ 
cal and spherical coordinates, differentiated only by sub¬ 
scripts would hopefully help in visualizing the complex 
geometry of skew ray tracing. 

In future papers, we shall extend our work on skew 
rays to ray tracing in circular, conical, and toroidal fibers, 
which may be step-index, multistep-index, or graded in¬ 
dex. 
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